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ABSTRACT 

We test the hypothesis that the temperature of the cosmic microwave background is consistent with 
a Gaussian random field defined on the celestial sphere, using de-biased internal linear combination 
(DILC) map produced from the 3-year WMAP data. We test the phases for spherical harmonic modes 
with ^ < 10 (which should be the cleanest) for their uniformity, randomness, and correlation with those 
of the foreground templates. The phases themselves are consistent with a uniform distribution, but 
not for £ < 5, and the differences between phases are not consistent with uniformity. For £ = 3 and 
£ — 6, the phases of the CMB maps cross-correlate with the foregrounds, suggestion the presence of 
residual contamination in the DILC map even on these large scales. We also use a one-dimensional 
Fourier representation to assemble aim into the AT^((^) for each £ mode, and test the positions of the 
resulting maxima and minima for consistency with uniformity randomness on the unit circle. The results 
show significant departures at the 0.5% level, with the one-dimensional peaks being concentrated around 
if = 180°. This strongly significant alignment with the Galactic meridian, together with the cross- 
correlation of DILC phases with the foreground maps, strongly suggests that even the lowest spherical 
harmonic modes in the map are significantly contaminated with foreground radiation. 

Subject headings: cosmology: cosmic microwave background — cosmology: observations — methods: 
data analysis 



1. INTRODUCTION 

' Since the release of the 1-year Wilkinson Microwave 
Anisotropy Probe {WMAP ) data (Bennett et al. 2003b,c; 
' Spergel et al. 2003; Hinshaw et al. 2003b; Komatsu et al. 
, 2003), great eflForts have been made to search and detect 

■ various possible forms of non-Gaussianity of the cosmic 
I microwave background (CMB) temperature fluctuations 

(Chiang et al. 2003; Coles et al. 2004; Park 2004; Eriksen 

■ et al. 2004c; Vielva et al. 2004; Copi, Huterer & Starkman 
; 2004; Cabella et al. 2004; Hansen et al. 2004; Mukherjee & 
'. Wang 2004; Larson & Wandeh 2004; McEwen et al. 2005; 
; Tojeiro et al. 2006; Dineen & Coles 2006). 

Three years of data are now available and methods 
• of foreground cleaning have also been improved, so the 
WMAP team have produced a new 'de-biased' version 
of their internal linear combination map (henceforth the 
DILC), which is claimed to be suitable for analysis over 
the full sky for spherical harmonic modes up to £ < 10 
(Spergel et al. 2006; Hinshaw et al. 2006). The statistics 
of these low multipole modes provide valuable information 
with cosmological significance, particularly concerning sta- 
tistical isotropy (or lack of it), demonstrated by the CMB. 

The 3-year data analyzed by the WMAP team is 
claimed to be Gaussian. Non-Gaussianity, if detected, 
could result from primordial origin (Bartolo et al. 2004), 
possible foreground residues left over after cleaning (Nasel- 
sky, Doroshkevich & Verkhodanov 2003, 2004; Dineen & 
Coles 2003; Naselsky & Novikov 2005), and/or correlated 
instrument noise. One will have to be very cautious if any 
non-Gaussianity is detected before attributing it to a pri- 
mordial origin. The concept of internal linear combination 
to obtain a reasonably clean CMB signal is to tune a set of 



weighting coefficients in order to minimize the variance of 
the foregrounds. Note that the variance is minimized but 
not eliminated. The possibility that foreground signals re- 
main in the DILC map to a significant extent is therefore 
something that should be carefully tested. 

In this paper we apply a series of stringent tests of the 
Gaussian hypothesis based on the behaviour of the phases 
of the spherical harmonic modes in the data. Since these 
phases are highly sensitive to the morphology (Chiang 
2001) of the temperature pattern, comparison between the 
phases of the DILC and the derived foregrounds should 
give a sensitive indication of the presence of contamina- 
tion. Phase information can also be used to diagnose 
departures from statistical homogeneity over the celestial 
sphere. 

Using spherical harmonic phases in statistical tests in- 
volves some subtleties. For one thing, they are not ro- 
tationally invariant. In other words a different choice of 
z-axis leads to a different assignment of phases for the 
spherical harmonic modes of the same pattern. This can 
be dealt with in a number of ways (Coles et al. 2004), 
but in the present context we choose to fix our coordinate 
frame as that which makes most sense given the probable 
behaviour of the foregrounds. In all the following we as- 
sume a Galactic coordinate system; all phase information 
is interpreted relative to this preferred frame. In doing this 
we attempt to ensure that the detection of non-uniformity 
or non-randomness in the phases can be interpreted more 
simply in terms of Galactic foregrounds. Issues such as the 
claimed north-south asymmetry (Eriksen et al. 2004b) and 
the alignment of multipoles (Tegmark, de Oliveira-Costa 
& Hamilton 2003; de Oliveira-Costa et al. 2004; Schwarz 
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et al. 2004; Land & Magueijo 2005), which seem to persist 
in the WMAP 3-year data, are also measured in Galactic 
coordinates, so we hope to shed some additional light on 
these peculiarities. 

Owing to the visual similarity of 1-ycar ILC and 3-year 
DILC maps, we will subject both maps to our analysis. 
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Fig. 1. — The for t < 10 are plotted on the Argand plane. 
The amplitudes are in unit of /xK and only m > modes are dis- 
played. The signs o, -|-, A, □ and * represent ^ = 1 — 5 on the left, 
£ = 6 — 10 on the right panel, respectively. Note that 12 out of 15 
phases of ^ < 5 are in the first two quadrants, among which six are 
clustered near or tt. 

2. THE GAUSSIAN RANDOM HYPOTHESIS AND THE CMB 

The statistical characterization of CMB temperature 
anisotropies on a sphere can be expressed as a sum over 
spherical harmonics: 



AT{e,(p) = ^ ^ aemYem{0,<p), 



(1) 



=0 m=-e 



where the Yf,m{0,ip) arc spherical harmonic fmictions, de- 
fined in terms of the Legendrc polynomials using 



Yem{e,'fi) = {-iy- 



I {2£ +!){£- my. 



Pem{cos 0) exp{im<p) 

(2) 

and the a^m are complex coefficients which can be ex- 
pressed with aem = \aem\ G^p{i<f>em) ■ In standard cosmo- 
logical models (i.e. those involving the simplest forms of 
inflation) these fluctuations constitute a realization of a 
statistically homogeneous and isotropic Gaussian stochas- 
tic process, or random field, defined over the celestial 
sphere (Bardeen et al. 1986; Bond & Efstathiou 1987). The 
formal definition of such a Gaussian random field requires 
that the real and imaginary parts of the aim are inde- 
pendent and identically distributed according to a Gaus- 
sian probability density, so that the moduli \aim\ have a 
Rayleigh distribution and the phases (j)e,n are uniformly 
random on the interval [0, 27r]. The Central Limit Theo- 
rem virtually guarantees that the superposition of a large 
number of harmonic modes will tend to a Gaussian as long 
as the phases are random, and this furnishes a weaker def- 
inition of Gaussianity. Moreover, statistical isotropy in 
general manifests itself in phase properties. Because of 
the importance of phases in both these definitions, we fo- 
cus on their measured properties as probes of departures 
from Gaussianity. Many methods have been proposed to 
test the random phase hypothesis (Chiang, Coles & Nasel- 
sky 2002; Chiang, Naselsky & Coles 2004; Dincen, Rocha 
and Coles 2005; Stannard & Coles 2005; Naselsky et al. 
2004; Chiang & Naselsky 2006). 



In Fig.l we plot on the Argand plane the agm of the 
DILC map for ^ < 10 (amplitudes |a^m| in unit of /xK.). 
Because of the conjugate properties of the aim for a real 
sky signal, we plot only a^m modes with m > 1 and omit 
all m = modes. Note the apparent non-uniformity of the 
phases for £ < 5. 



3. TESTING THE RANDOM PHASE HYPOTHESIS 

We use Kuiper's statistic (Kuiper 1960) (KS) to test 
on the random phase hypothesis. The KS can be viewed 

as a variant of the standard Kolmogorov-Smirnov test, 
designed to cope with circular data. The standard 
Kolmogorov-Smirnov statistic is taken as the maximum 
distance of the cumulative probability distribution against 
the theoretical one: D = max_oo<2,<oo \Sn{x) — P{x)\. 
For a circular function, however, one needs to take into 
account the maximum distance both above and below 
the theoretical probability P{x). Accordingly, we define: 
D± = max_oo<a;<oo ±('S'Ar(a;) — P{x)). Then the test 
statistic is 



V = D++D_. 



(3) 



In standard frequentist fashion, we define the significance 

level a (or p- value or "size") for our "mfll" hypothe- 
sis (of uniform randomness of the phase angles) as the 
probability of the measured value of V arising under the 
null hypothesis. In this case, a can be calculated from 

a = QKuiper(V^[V^ + 0.155-F0.24/ViV]) and QKuiper(A) = 

2 (4j^A2 - l)e-2-''^' , where N is the number of data 
points. 

We perform three statistical tests based on this general 
approach: on the uniformity of phases (i.e. consistency 
with a uniform distribution on the interval [0, 27r]); on the 
randomness (i.e. independence) of phases by taking dif- 
ference of phases with fixed separation {A£, Am); and on 
the cross-correlation of each £ between DILC and the fore- 
grounds by A(j)f^ = (j)"^^ — (ff^. In each case the resulting 
angles should be random: the difference between any two 
random angles is itself a random angle. 





4>em of ^ < 5 


(pern of 6 < ^ < 10 


all (t>em, of ^ < 10 


ILC 


38.74% 


90.91% 


73.55% 


DILC 


30.63% 


98.09% 


66.14% 



Table 1 

The significance level a in accord with the uniform 

DISTRIBUTION HYPOTHESIS OF THE PHASES OF THE 1-YEAR 

ILC and 3-year DILC maps. 

We first test the uniformity of the phases for two groups 
{1 < £ < 5), and {Q < £ < 10) and for all phases claimed 
to describe clean modes (i.e. i < 10). The phases of the 
TO = modes are excluded in all our test. Table 1 shows 
that for the 3-year DILC, while the phases are consistent 
with a uniform distribution (a ~ 0.98) for 6 < £ < 10, and 
for ^ < 5 (a ~ 0.31). This is consistent with the behaviour 
of the phases seen in Fig.l: the phases seem to be concen- 
trated in the first two quadrants. This is however, only 
significant at the 31% level, which is barely la. Overall, 
therefore, the phases of the 3-year DILC map for £ < 10 
are consistent with uniformity {a = 0.66). 
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Fig. 2. — Significance levels of randomness with separation 
(A£, Am). The thick and thin lines are from the 3-year DILC and 
1-year ILC maps, respectively. 

Now we check the randomness of the phases is tested 
by the defining a set of difference A(/)'^'''^(A£, Am) = 
^f+A£,m+i\m ~ <^^m ■ Fig.2 wc show the significance 
levels for randomness between phases with different I and 
m. There are three separations that show significant de- 
partures from uniformity (a = 0.047, 0.0421 and 0.0053 
for (A^,Ato) = (0,1), (0,3) and (1,2), respectively. This 
corresponds to significant coupling of the phases. It is 
complicated to understand coupling across both £ and m 
for (A£, Am) = (1, 2) because of the lack of rotational in- 
variance described earlier. Nevertheless, because the first 
two examples involve coupling between azimuthal numbers 
m within each we plot in Fig. 3 the sequences of phase 
difference A0'^'''^(0, 1) and A0'^'''^(O, 3) on a unit circle (or 
exp(iA0'^'''^) on an Argand plane). For both distributions 
one can see the deficits around A(/) = that cause the 
significance levels for the cases to appear below 5%. We 
also include (0, 2) in Fig. 3 for comparison, which has the 
apparent tendency of phase differences to avoid tt, though 
this is not significant at the 5% level. 

We would like to emphasize that the statistics derived 
from these 35 separations shown in Fig. 2 should not be 
treated as a statistical ensemble. For example, primordial 
cosmological magnetic field induces and supports vortic- 
ity or Alfven waves, which induce in CMB anisotropics 
with correlation between ai+i^m and a£_i i.e. A£ = 2 
(Durrer, Kahniashvih & Yates 1998; Chen et al. 2004). 
Therefore, examining A£ = 2 correlation alone is qualified 
as an independent method. Another example is that sym- 
metric signals defined on the Galactic coordinate system 
with respect to the Meridian [tp = 0) on a sphere induce 
correlations between A£ = 4 (Naselsky & Novikov 2005; 
Naselsky, Novikov & Chiang 2006). Therefore, each of the 
(A£, Am) should be treated as a separate non-Gaussianity 
test like bispectrum and trispectrum. . . etc.. 

4. CROSS-CORRELATION BETWEEN THE DILC AND 
FOREGROUND MAPS 

Since the DILC map is obtained from an internal com- 
bination of the frequency maps, some foreground residues 
might be left unsubtracted. Based on the assumption that 



CMB signal should not correlate with the foregrounds, and 
that the characteristic of phases refiect the morphology 
of the CMB anisotropy pattern, we also test the cross- 
correlation of phases for the DILC and the WMAP 3-year 
foreground maps at K, Ka, Q, V and W channels. The 
foreground maps we test are the sum of the synchrotron, 
free- free and dust templates. We take the phase difference 
^^fm = ^fm ~ ^Jm foi' ^ach t, assumiug such A(/)£,„ at 
each i should be uniformly distributed. In Fig. 4 one can 
see that for £ = 3 and 6 the DILC phases have correlation 
with the foregrounds with significance around 10%. One 
particular point about the quadrupole is easily seen from 
Fig.l: two of the three quadrupole phases are near and 
TT. Note also that, for £ = 6, in Fig. 3 (shown in red dia- 
mond sign), the phase differences for Am — 1 are strongly 
clustered. 
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Fig. 3.— Phase diff'erence for Am = 1,2 and 3 for all <l>f^, i < 10, 
1 < m < ^ plotted on a unit circle. The significance levels of depar- 
tures from uniformity are 4.70%, 20.85% and 4.21%, respectively. 
The signs o, -|-, A, □ and * in black represent phase difference within 
^ = 1 — 5, respectively, and in color within £ = 6 — 10, respectively. 



5. EXTREMA STATISTICS 

The Am phase coupling in each £ leads to a departure 
from Gaussianity in the resulting signal. Following Chi- 
ang & Naselsky (2006) we represent the effect of phase 
coupling by assembling the aim (now only single variable 
m) with an inverse Fourier transform: 



(4) 



where the negative m in the sum are replaced with a^^^ to 
ensure that ATg is real. The morphology of the signal ob- 
tained by this method is similar to the signal AT{9, ip) = 
X]m=-£ '^eniYim^O , if) obtained for each £ by summing over 
all 9 onto ip axis on the map. One should note the following 
subtleties of such comparison. The map can be written 



ATf 



I2£+1 



aioPioicose) 



m— 1 



'{£-my. 
{£ + m)l 



aijn\cOs{m(p + (j)ira)Peni{cOS0) 



■(5) 



The integration of the associated Legendre function is 
zero for odd £ + m. When one compares the ID curve 
with the composite map for odd £, only the odd m 
contribute to the integration while the (—1)™ reverses 
the signal as seen from simple inverse Fourier transform 
\aem\ cos{mp> + (pt 

Such ID construction obviously loses information in one 
of the available dimensions, but despite the fact that the 
a^m are spherical harmonic coefficients, the statistics reg- 
istered in the one-dimensional complex aim should still 
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manifest themselves in the ID AT curves we create. In 
Fig. 5 we plot AT{(p) summing from the DILC aem (top), 
and from the whitened DILC: aim/\aim\ (bottom). If the 
signal is Gaussian, the locations of the highest and lowest 
peaks should randomly distributed in ip between —180° 
and 180°. We plot in Fig. 6 the distribution of these ex- 
trema locations in a unit circle. 
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Fig. 4. — Significance levels for the cross-correlation of phases at 
each £ between the DILC (thick lines) and the WMAP foreground 
maps, from top to bottom, K, Ka, Q, V and W channels, respec- 
tively. For comparison, the thin lines are those between 1-year ILC 
and 1-year foreground maps. 

In Table 2 we list the significance of the distribution of 
the extrema locations. Not only do they show low signif- 
icance levels of random distribution in (p, the peak loca- 
tions also cluster in ip = 180°. In Fig. 4 the i = 3 and 
6 show significant cross correlation with the foregrounds. 
Wc therefore test the peak distribution by excluding the 
four extrema belonging to these two modes. Although p- 
values increase slightly, the results are still significant at a 
level l)d()w 5%. 
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0.552 % 


excluding peaks of £ = 3, 6 
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Table 2 

sicjnificance levels of the uniformity of the 
distribution of the extrema locations in tf. 
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Fig. 5. — The AT^((^) distribution for each £ from the DILC 
map by assembling J],^ a^^'^exp(im<p) (top panel, in which the 
unit of y axis is /iK) and by assembling only the phase part: 
Sm '^-'^P(*''^fm'^) '^-'^P(*™''') (bottom). In order to match the conven- 
tion of Galactic longitude coordinates (for comparison with, e.g. 
Fig. 14 in Hinshaw et al. (2006)), the ip axis is plotted reversely. 
Note that on comparison, due to the properties of spherical har- 
monics, peaks on the maps for odd I correspond to troughs on the 
ID curves. In each figure we indicate the locations of the extrema. 

6. CONCLUSION 

In this letter we test the Gaussian random hypothesis 
of the CMB temperature anisotropics. The behaviour of 
the 3-year DILC does not differ strongly from that of the 
1-year ILC version: all the famous peculiarities still ex- 
ist, as mentioned in Spcrgcl ct al. (2006). Wc have found 
in this paper that this also extends to the behaviour of 
the spherical harmonic phases. In particular, we find that 
phase differences (which should be uniformly distributed), 
tend to avoid the region of the complex plane close to 
/S.(f) = 0. We also find that the phases for £ = 3 and £ ^ 6 
are significantly correlated with those of the foreground 
maps. We also test the real space alignment of the re- 
sulting features using the temperature extrema resulting 
from a Fourier summation for each £. The resulting peaks 
are indeed concentrated at opposite to the center of our 
Galaxy, i.e. I = 180° of Galactic coordinates, with respect 
to which the phases are themselves defined. 
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On the basis of these results we reject the null hypothesis 
that the modes with I < 10 are a realization of a statisti- 
cally homogeneous Gaussian random field at a significance 
level better than 5%. We also infer that the origin of the 
observed departures is consistent with being some form of 
Galactic foreground. 

Of course it is possible the apparent alignment between 
CMB temperature pattern and galactic foreground mor- 
phology is simply fortuitous. It could be that large-scale 
anisotropics in both line up accidentally. If this is the 
case then we just happen to live at a place in the Universe 
where our past light cone presents us with this coincidence 
and we will just have to cope with it; all future diagnostics 
of foreground contamination will have to incorporate this 
alignment as conditioning information. We believe how- 
ever that it is important to take such coincidences very 
seriously until we know for certain that is all they are. 
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Fig. 6. — The distribution of the extrcma locations of AT£(ip) as 
shown in Fig. 5. The angle of the unit circle denotes the location 
if. The left panel is for the DILC ATi{ip), and the right for the 
whitened DILC AT((ip). The + sign denotes maxima and o sign 
minima. 
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